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Abstract 

We calculate the one loop gluon vacuum energy in the background of a 
color magnetic vortex for SU(2) and SU(3). We use zeta functional regu¬ 
larization to obtain analytic expressions suitable for numerical treatment. 
The momentum integration is turned to the imaginary axis and fast con¬ 
verging sums/integrals are obtained. We investigate numerically a number 
of profiles of the background. In each case the vacuum energy turns out 
to be positive increasing in this way the complete energy and making the 
vortex configuration less stable. In this problem bound states (tachyonic 
modes) are present for all investigated profiles making them intrinsically 
unstable. 


1 Introduction 

Color magnetic vortices or center-of-group vortices play at present an important 
role in the discussion of possible mechanisms for confinement in QCD Jl|, [2. d||. 
There are recent investigations of vortex like configurations in QED [j|] and in 
QCD ||. Given the actual interest in this topic we consider in the present paper 
the case of Abelian color magnetic vortices with arbitrary profile for SU(2) and 
SU(3) gluodynamics. Following the technique developed in 0 for a magnetic flux 
tube and earlier in [|7| for a spherically symmetric background field we write down 
analytic expressions for the vacuum energy well suited for numerical investigation. 
Using zeta functional regularization we express the vacuum energy in terms of 
the Jost function of the related scattering problem, and by separating a piece of 
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the asymptotic expansion of the Jost function we obtain fast converging sums 
and integrals. 

At present, several calculation schemes are known. The oldest and perhaps 
most intuitive one is to sum over phase shifts or, what is equivalent, over mode 
densities. For recent applications see 0 and [j|]. This method is numerically 
demanding because of the oscillatory behavior of the phase shifts as function of 
the radial momentum. Another method that had been proposed in § uses the 
Jost function and a momentum integration turned into the imaginary axis in this 
way avoiding oscillatory integrands (for a detailed description of these methods 
see also 0). A completely different approach is that of [110, O using world line 


methods. Its main advantage is that it does not require separation of variables. 
Recently, this method was applied to a magnetic string and numerical results 
similar to those in [jBJ have been obtained (for a comparizon see [12 . 


The issue of renormalization of vacuum energies has been much discussed in 
the past and different schemes are in use. Finally, all differences are settled by 
an appropriate normalization condition and the many schemes are different ways 
to the same destination. However, this is not always easy to see. Let us first 
consider theories with a massive quantum held, a spinor for instance. Here the 
mass of the held provides the opportunity to formulate the ’large mass’ normal¬ 
ization condition stating the vacuum energy must vanish if this mass becomes 
large [|1^, [14|. This condition was initially formulated in connection with the 
Casimir effect for massive helds |T5|. 


In |T6| 


it was shown to be equivalent to 
the no tadpole normalization condition known in quantum held theory. This 
condition is directly related to the heat kernel expansion which is equivalent 
to the expansion of Greens functions in inverse powers of the mass |T7]]. It is 
also evident that this condition is in line with the charge renormalization ap¬ 
pearing from vacuum energy calculations in QED. Another picture we observe 
with massless quantum helds. There is no such normalization condition causing, 
for instance, the known problems with the vacuum energy of a dielectric sphere 
[ |TJ|j . However, in the present paper we are concerned with QCD and the vac¬ 
uum fluctuation of the gluons. Although they are massless, in QCD there is the 
commonly accepted coupling constant renormalization involving the phenomeno¬ 
logical quantity A qcd- We will use this scheme. The only technical problem is 
that this scheme is frequently used together with the Pauli-Villars regularization 
while we use the zeta functional one. But the relation between both schemes can 
easily be written down, see Sect. 2. 

The paper is organized as follows. In the next section we collect the general 
formulas for the color vortex and its vacuum energy. Following 0 we write down 
the condition for the minimum of the complete energy. In the two subsequent 
sections we consider in detail the finite and the asymptotic parts of the vacuum 
energy for the given configurations. In Sec. 5 the numerical results are repre¬ 
sented. Sec. 6 contains a discussion and some technical details are given in the 
Appendix. 
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We hope to have represented all formulas in such detail, or given the corre¬ 
sponding references, that the reader is able to repeat the calculations. 

2 The color magnetic vortex and its vacuum en¬ 
ergy 

We consider a non Abelian gauge field with standard action 

S = 4? / ^ (W) 2 ’ « 

The corresponding potential, A“(x), is expanded around a background, A^ bg ^(a;), 

a;(x) = a^(x) + a;^(x), ( 2 ) 

where A^ qu \x) are the quantum fluctuations. The background is a center-of- 
group vortex, i.e., it is an Abelian static and cylindrically symmetric configura¬ 
tion. The corresponding color magnetic field is parallel to the third axis in both, 
color and space, 

B a (x ) = (W 3 e z (3) 

r 

where /i{r) is a profile function, e z is the unit vector parallel to the z-axis and 
5 a ’ 3 is the Kronecker symbol. The classical energy of this field meant as energy 
density per unit length with respect to the z-axis is 



The profile function will be chosen vanish at r = 0. Its first derivative must 
vanish too otherwise the classical energy is not finite. For r —> oo it must tend 
to a constant which determines the flux. We consider examples with //(oo) = 1. 

The fluctuation fields A“( qu )(a;) are quantized in a standard way. We use the 
background gauge condition as in || and have gluon and ghost contributions. 
The vacuum energy due to these fields is 

c _ c gluon , c ghost 

£ quant ^ quant ' ^quant- 

In general, a vacuum energy is given by 

f .-t-VA 24 

°quant ~ az 2 Z^ C (n) > 

(n) 

where e( n ) are the corresponding one particle energies and the sum over (n) is 
over all quantum numbers present for the given held. The parameter s is the 


(5) 

( 6 ) 
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regularization parameter of the intermediate zeta functional regularization which 
we use. We have to take s > | in the beginning and to put s = 0 in the 
end. We ignore the arbitrary parameter /i which is usually introduced in this 
scheme. The sign is +1 for gluons and -1 for ghosts. For both types of particles 
there is a sum over color and polarization. For the ghosts there is an additional 
factor of 2 because there are two fields, c and c. In the considered background, 
variables separate in the equations of motion for the fluctuating fields. One can 
separate the dependence on the z-coordinate and the angular dependence in the 
plane perpendicular to the z-axis introducing p z and the angular momentum v 
(u = 0, ±1, ±2,...) accordingly. After that the one particle energies read 


e W = 



(7) 


where AyO' are the eigenvalues of the residual radial equation, 


Id d (y — a c p{r)) 2 

ry » |^/^i Qj, ry* 2 


2 a c f3xp\r) \ 
r ) 




( 8 ) 


Here we assume the radial quantum number n to be discrete, see below. This 
equation emerges after diagonalization in color space. The corresponding color 
eigenvalues are 


OL c 


{0, -1,1} for SU(2) 

{0,0,-|,—A, |, A,—1,1} forSU(3). 


(9) 


For the ghost field we have to put f3\ = 0 in Eq.(||). For the gluon the polarizations 
/3\ are 

= { 0 , 0 ,- 1 , 1 }. ( 10 ) 

In the vacuum energy, the two gluon polarizations with f3\ = 0 cancel the ghost 
contribution. So we are left with the two gluon polarizations /3\ = ±1 and no 
ghosts. 

In the next step we integrate over p z by means of 


f°° dpz 

J —oo 27T 


(p^ + A 2 )* ' = A 2 < 1 -‘> 


r<> - 1 ) 

2-/HT(s - i; 


= ^2(1-.) C, 


47TS 


with C s = 1 + s(2 In 2 — 1) + ... and arrive at 


00 0/1 \ 

£,u. n « = fEEEW (ii) 

Z c,A v=-co n 

With this representation and equation (p|) the whole problem is reduced to a 
number of scalar ones. Starting from here, and for pure technical convenience we 
introduce an auxiliary mass m by means of 

Aft* 1 - (A<y )2 +m 2 )' /2 (12) 
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with m —> 0 before the end of calculations. 

It remains to transform the discrete sum over the radial quantum number 
n into an integral. As the initial problem is defined on the unbounded plane 
perpendicular to the z-axis this quantum number is continuous. We make it 
discrete by first imposing a large ’box’, i.e., a cylinder at some large radius R 
with some boundary conditions. We follow || where this procedure is described 
in detail, but we discuss all subtle specifics of the given case. So, omitting all 
indices, let (j){r) be the regular solution of the radial equation on the whole axis. 
It is defined as that solution which at r —> 0 turns into the solution for /i{r) = 0 
which is simply the Bessel function J u (kr). For large r, <p(r) has the asymptotic 
expansion 

</>{r) ~ \ (fu(k) H^\kr) +J v (k) H^\kr)\ . (13) 

r—>oc Z X 7 

The coefficients f u (k ) and its complex conjugate are the Jost functions. The Han- 
kel functions are also solutions of the radial equation (|8|) for /i(r) = 0. However, 
we have to take into account air) ~ 1 and should use the Hankel functions 

' v 7 r—KXD 

H^Za c (kr) instead. For the definition of the Jost function, this doesn’t matter 
because the leading term in the asymptotic expansion of the Hankel functions for 
large argument doesn’t depend on the index. 

Following ||, for large /?,, the boundary condition <f>(R) = 0 defines the dis¬ 
crete eigenvalues k —> X n . This allows to write the sum over n as a contour 
integral encircling the spectrum with the logarithmic derivative of 4>{R) generat¬ 
ing the necessary poles. After that we tend r to infinity and drop a contribution 
which does not depend on the background (the so called Minkowski space con¬ 
tribution). In the heat kernel terminology this is to drop the coefficient cio’s 
contribution. Finally turning the integration contour toward the imaginary axis 
we arrive at the representation 


E K + 


m 


i—s sm ns 


—n 


r\ 

dk (k 2 — m 2 ) —In fJik). (14) 
v 'ok 


Putting these things together, we get for the vacuum energy (|TT|), 


= ^51U> 


(15) 


where now f^ c,x \ik) is the Jost function corresponding to Eq.(^) and we kept 
only first order corrections for s —> 0. This is a very convenient representation 
for the vacuum energy. The main advantage is that the integration goes over the 
imaginary axis avoiding oscillating functions. In addition we remark that bound 
states, if present in the spectrum, are included in this formula automatically!;. 
Bound states manifest themselves as zeros of the Jost function on the positive 

1 This is discussed in detail in |l8j and used in numerical examples in Q 
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imaginary axis. If the binding energy K t is smaller that the mass m of the fluctu¬ 
ating field the corresponding poles of the logarithmic derivative are outside the 
integration region in Eq. flT5|) and the vacuum energy is real. In opposite case, 
for m < Kj, the integration runs over the pole (it must be bypassed on the right) 
and an imaginary part appears. Using / = in we obtain 

3£quant = ~ > ( 16 ) 

° c, \ V i 


where the sums go over those values of the quantum numbers for which a bound 
state is present. For instance, i is limited by m < In our case we have rn — 0 
and any bound state contributes to the imaginary part of the vacuum energy. 
We note that the imaginary part does not contain ultraviolet divergences and so 
that we could put s = 0 here. The real part of the vacuum energy is then given 
by Eq. (]T5|) taking the integration over k as principal value. 

The vacuum energy given by Eq. fll5D still contains known ultraviolet divergen¬ 
cies. Their structure can be best discussed in terms of the heat kernel coefficients. 
The contribution of the coefficient a 0 we already dropped. Due to the gauge in¬ 
variance the coefficent Gq is zero. The next coefficient is 02 which is here given 
by 


-22 N 


dr 


ci 2 — 


-7T 


Arf 


(17) 


At once it is the last coefficient which contributes to the ultraviolet divergencies. 
By writing its contribution to the vacuum energy separately (for a derivation of 
the corresponding formula see J7|) we define by means of 


S, 


-Cl 2 


quant 


32tt 2 


2 - 2 In 


m 


4 - P ren 
1 ‘“'quant 


(18) 


the renormalized vacuum energy £™( nt . bi general, after separating the ultraviolet 
divergence, the remaining finite (in the limit of removing the regularization) part 
is not uniquely defined. In the form given here it contains contributions from 
the heat kernel coefficients a n with numbers larger than 2 only. Because the heat 
kernel expansion provides an expansion in inverse powers, m 4 ~ 2n , of the mass of 
the fluctuating held, £ q ™ nt goes to zero for m —► oo. In problems with a massive 
fluctuating held we used this property as normalization condition defining in this 
way a unique renormalized vacuum energy see also [2I)j]. For massless 

helds this does not work. In the case considered in this paper the underlying 
theory is QCD. Here the renormalization is done in the Pauli-Villars scheme by 
a redefinition of the bare coupling constant in Eq.([l|) according to 


1 _ 111V M 
g 2 24tt 2 11 Aqcd’ 


(19) 


where M is the Pauli-Villars mass and N is the order of the gauge group, SU(N). 
We denote for a moment the regularized vacuum energy by ^ qU ant(s, rn), indicating 
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the dependence on the regularization parameter s and on the mass m explicitly. 
In the Pauli-Villars scheme the vacuum energy is given by 


c pv 

^ quant 


^quant (5, Tfl) <?quant (5, M ) 


( 20 ) 


for m —► 0 and M —> oo. Using the representation of the vacuum energy given 
by Eq.(|T8|) as a result of the heat kernel decomposition we obtain 


c PV 

^ quant 


~a 2 /l 
32tt 2 Vs 
-02 / 
32vr 2 V 


777 \ 

-2-2 In —j + £^ nt (s,m) 

1-2-2111 y) - f“„t(M'/). 


( 21 ) 


Here a number of contributions proportional to a 2 cancel, the pole term for 
instance^, and in the limit we are left with 


<?pv 

^quant 


— o 2 
167T 2 


111 M + ^2 111 m + ^quant M ■ 


87T 2 


( 22 ) 


Here we have taken into account that £^ nt (s, M) vanishes for M —> 00 . Next, 
the contribution of In M is canceled from the corresponding contribution from 
the classical energy by means of ([19|). Finally, in ([22jj) the limit m —> 0 gives a 
finite result. This follows because £ quan t as given by Eq. ([18|) is finite for a massless 
theory. 

For the total energy we obtain 


£ 


£ 


class 


+ £ PV 

1 ^quant 


11 IV 

247 


In Aq CD 



p'(r) 


2 1 crenPV 
' ^ quant 


with the notation 


crenPV 
^ quant 



lnm + ^ant 



\m=0 


(23) 


(24) 


for the vacuum energy renormalized in the Pauli-Villars scheme. 

Eqs. (|23|) and (^) are the final answer for the complete energy of the vortex 
in QCD. It will be calculated for a number of profile functions in the subsequent 
sections. Because the vortex itself is unstable as its classical energy is positive a 
physical interest is in configurations lowering the complete energy. To find those 
a dimensional analysis can be made following ||. First we note that the profile 
n(r) is a dimensionless function of a variable with the dimension of a length. 
Therefore it must contain an intrinsic scale r 0 so that in fact we have to use 


P 


r 

p ( — 

r 0 


2 Starting from here we can put s = 0 in fqJVnt- 


( 25 ) 
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Now we introduce the notation 

POO rjrp 

A = ir — //(r) 2 (26) 

Jo r 

to be used in Eq.(|23|)- Note that A is dimensionless and a factor of 1/rg had 
been taken out of the corresponding integral in Eq. (|23|) . Next we note that 
£quant V > Eq. (|4]) has for a background profile as given by Eq.(|^) the following 
dependence on r 0 , 


crenPV 
^ quant 


11 N 

247r 2 rg 


A In r 0 + -$£, 


renPV 


2‘-'quant | ro= l’ 


This is for dimensional reasons. We introduce by means of 


(27) 


111 D 


247T 2 prenPV 

11NA C ' quant ko=i 


(28) 


the notation B , which allows to represent the complete energy in the form 

„ UN A , . . 

£ ~ 11 ( r ° A Q CDjB ) • ( 29 ) 

After this rewriting the dependence of the complete energy on the scale ro is 
present in the shown places only. It is easy to find the minimum as a function of 
rg. After trivial calculation the minimum appears at 


1 y/e 
r ” 1 "" “ Aqcd IP 


(30) 


where e is the basis of the natural logarithms. The energy in the minimum is 


£ . = -A 2 -- -. (31) 

Imm QCD 247T 2 2e 

This is the quantity to be compared for different profiles. Here the sign of the 
vacuum energy enters as follows. Primarily, the sign of the complete energy in the 
minimum, £| , is always negative. This is a consequence of the sign of 02 resp. 

of the first coefficient of the beta function. Secondly, the sign of in Eq.([27|) 

is dominated by In r 0 in the first contribution on the rhs of Eq. (^) and depends 
in this way on the scale. The second contribution, Tq™(( t V | r() _|, ma Y be positive 
or negative which will be a result of the calculations for a specific profile p(r). 
If the sign of t is positive (negative) we have B < 1 (B > 1) and this 

contribution increases (decreases) the complete energy. At once, B influences the 
core radius in the minimum, r 0 | min , making it larger (smaller) for B > 1 (B < 1). 
We will see below in all considered examples we have B > 1. After this discussion 
of the scale dependence in all following calculations (Sec. 3 and 4) we put r 0 = 1. 
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3 The finite part of the vacuum energy 


The regularized vacuum energy is given by Eq. (|15|) and we need to represent it 
in the form of Eq.(P^j), i.e we need to separate the pole part in s. Eq.(|T5|) is well 
defined for s > § and we need to perform the analytic continuation to s = 0. In 
order to do that we add and subtract a piece of the uniform asymptotic expansion 
of the logarithm of the Jost function which we define by 




ln/,(zfc)~ln/r(zfc) + 0(^ 

h) 


(32) 

for v, k 

—> oo with z e 

e * fixed. In 

V 

this way we represent 





? — 
^ quant 

£ { + £ as + £l + £ 0 as 


(33) 

with 







£ { = 

-I OO 

r oo / 

1 dk (k 2 - 

m ' 

ro2 ) It ( ln/ ” c,A,( 

:«) - in /y- A) 

“(«)). 

(34) 








and 








£ as = Cs 
—8n 

v roo 

EE/ 

e,A Jm 

dk [k 2 — m 

— In f( c > A )as/ 

dk Ju [ 

ik ), 

(35) 


i /#0 


where we in the Jost functions again indicated the dependence on color and on 
polarisation. In Eq.(^) we have put s = 0 because there the sum and the 
integral are convergent. This is ensured by the choice of the asymptotic part of 
the Jost function. It can be checked that both, the integral and the sum are 
made convergent by subtracting the uniform asymptotic expansion. This fact is 
non trivial and doesn’t hold for any background potential. For instance, it can be 
shown that for a /i(r) with a non vanishing first derivative in r = 0 the integral 
is not convergent. But in that case the classical part of the energy, Eq.(f|), is 
infinite. This we do not consider. 

In Eqs. (|34]) and fl35|) we left out the contribution from the orbital momentum 
v = 0 because it needs a separate treatment. Here we define the asymptotic part 
of the Jost function as that for large argument k , 


h(ik) ■- + O (j^J 


(36) 


for k 


oo. 


Accordingly we have 


and 



Al„ (/f’ A) («) -/h A, “(«)) , 

-m 2 ) 1- * A in/Ml-jjj,). 


(37) 


(38) 


9 





Again, in Eq.(|37]) we could set s = 0 because the integral over k is convergent. 

The asymptotic expanson of the Jost function can be obtained in the same 
way as in Ref. [[?!] and in Ref. IJ from iterating the Lippmann-Schwinger equation. 
While in 0 the background was scalar only two iterations were needed while 
in [16] for a vector background four had to be included. In the present case the 
background is a vector, hence we need four iterations. In |J the quantum field 
was a spinor one and here we have a vector held. So we need to modify the 
formulas of [6] to some extend. The basic formula is the Lippmann-Schwinger 
equation which appears if rewriting the wave equation (||) as an integral equation 
(we drop the color and polarization indices for a moment), 

r , , 

0(r) = 0°(r) + Yif (Mkr)H^\kr') - J v {kr')H^\kr)) P(r / )0(r / ) (39) 

o 

with the ’perturbation operator’ 

P{r ) = —2 ct c p(r) + a 2 ^{r) 2 — 2a c (5xrii'{r) (40) 


and the choice 0°(r) = J u [kr ) for ’free’ solution which just makes 0(r) defined 
by Eg. (fiO|) the regular solution. From considering this equation for r —> oo and 
comparing with Eq.([T3|) the following representation for the Jost function can be 
obtained, 

7T 7 dr 

Uk ) = 1 + yJ -HV(kr)P(r)<t>(r). ( 41 ) 

0 

The next steps go exactly as in Ref.||. The logarithm of the Jost func¬ 
tion which is the quantity we are interested in can be represented as given by 
Eqs.(32)-(35) in || with obvious substitutions. In these formulas one has to turn 
to imaginary k. After that the modified Bessel functions entering have to be sub¬ 
stituted by their well known uniform asymptotic expansion and the integrations 
have to be performed by the saddle point method keeping all the contributions up 
to the necessary order. This is a quite lengthy calculation which is best performed 
by a computer algebra package to avoid errors. The result is a representation of 
the asymptotic expansion of the logarithm of the Jost function of the form 

3 3n 7 dr t j 

ta/r («0 = ££/-w j - ( 42 ) 

n=lj=n J Q ' V 


with 


t = 


1 

\Jl + ( zr ) 2 


(43) 


(z = k/v). The coefficients X n j are local functions of p(r) and its derivatives. 
They are listed in the Appendix. It turns out that this is not the most convenient 
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form. It is possible to integrate in Eq.([l2|) by parts so that all contributions 
become proportional at last to one derivative of n(r). In particular this speeds 
up the convergence of the integrals at r —> oo which later on must be taken 
numerically. After integrating by parts the formula looks essentially the same 
except that the coefficients are different. We denote them by Z n j so that we 
obtain In fjf(ik) in the form 



where the new coefficients are listed in the Appendix too. 

For zero orbital momentum, v = 0, we need the asymptotic expansion of the 
Jost function for large argument. It can be obtained from Eqs.(32) and (33) in 
Ref. 0 by inserting there the asymptotic expansion of the Bessel functions for 
large argument. This is much easier than in the case of the uniform asymptotic 
expansion. This allows us to present these formulas here in order to illustrate the 
method. From the mentioned equations in || we have 

°r d r 

In / (1) = —Ukr)KJkr)P(r ) (45) 

J r 
o 


and 

°r dr °r dr' 

In = / — K%(kr)P(r ) / — I^(kr')P(r'). (46) 

Jr Jr' 

o o 

In Eq.([|5|) we use the expansion 


I v (kr)K v {kr) 


1 1 
2 kr 16 (kr) 3 ^ 


(47) 


In Eq. (|46|) it is sufficient to use the leading order from the expansion of the Bessel 
functions, K'l(kr) = exp(— 2kr) / (2kr) + ... and J 2 (Ax) = exp(2/cr)/ (2kr) + ... 
and to take the integral over r J by the saddle point method to leading order. The 
result is 

lll/,2, = -8F Jv Pir)2 - (48) 

0 

Taking both contributions together we define 


in /„“(«) = / 

o 



P(r) 
2 kr 


+ 


1 

( fc 2 + 1) 3/2 




(49) 


Here we substituted 1 /k 3 by 1 /{k 2 + 1) 3 / 2 . This does not change the asymptotics 
to the given order but avoids problems with the integration over k at the origin 
in Eq.(|37|) when we put m — 0 later on. 
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Now, in order to calculate the finite parts of the vacuum energy we need 
a method to calculate the Jost function numerically. Here we follow the ideas 
of Ref.JIiJ with necessary modifications. The Jost function f u (k) is defined by 
Eq. ( p~3|) . We start from defining a new function f u (k,r ) by means of 

0(r) = \ (f v (k, r) Hl 2 } ac (kr) + J u (k, r ) H^krj) . (50) 

From the definition, Eq. fll3|) , it is clear that 


U(k,r ) ~ fv(k). 

r—»oo 


(51) 


The speed of convergence depends on how fast the potential fi{r) approaches its 
value at infinity. It can be shown that the derivative of f u (k,r) with respect to r 
is proportional to /i(r) — 1 . 

For a numerical calculation of the wave function (f>(r) we can use the radial 
equation, Eq. (^|) and solve it with some Runge-Kutta procedure. For this end we 
need initial conditions at r = 0. Here it is necessary to introduce a new function 
by means of 

/ kr\ u <p(r) 

\YJ r(v+iy 



Given <p(r) is the regular solution it behaves for r —> 0 like the Bessel function 
J u (kr ) and the introduced function 99 (r) is finite at the origin, 


<P(0) = 1- (53) 

As the second initial value we need the derivative of ip(r) in r = 0. Because fi{r) 
vanishes there together with its first derivative, 

tp\ 0) = 0 (54) 


holds. 

In the numerical calculations we used Mathematica and the build in function 
NDSolve. Here it is useful to turn to a system of first order equation by means of 

•0V) = -g^±ly, (r) + + _ k2 \ 

r y r r 2 r j 

Also it is necessary to move the starting point of the procedure away from r = 0 
because otherwise NDSolve tries to divide <£>'(0) by r = 0. The initial data shifted 
by a small amount A can be obtained from a decomposition of the equations in 
powers of r and in the lowest nontrivial order we get 

V?(A) = 1, 

^(A) = A /i 2 (55) 
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with ji -2 = n"(r )| r=0 . Now it remains to express the Jost function, f v (k,r ), from 
Eq. (|50| ) in terms of the functions (p(r) and ^(r). In addition we need the Jost 
function for imaginary argument. So we have to substitute k —► ik in the above 
equations. An easy calculation gives 

f„(ik,r) = AA— ^rip(r)K y _ a (kr) + ip(r) (kK„ + ,_ a (kr) + ?-K„_ a (kr)^ , 

(56) 

where K v _ a (kr) are modified Bessel functions. In the process of solving the 
equations using NDSolve the true Jost function is approached quite fast. In all 
considered examples it was sufficient to take r = 10. 

By means of Eqs. (^) and (£J3|) we have all we need for a numerical calculation 
of £ { . In the actual calculations we used Eq. fl34|) in a modified form. We inte¬ 
grated by parts and change the integration variable from k to z = k/v. Finally, 
we put m — 0 because here this auxiliary mass is no longer needed. In this way 
the formula reads 

1 OO n QQ 

^ f = rEEM dzz (in fl c ’ X \iiyz) - In f^ as (ivz)) . (57) 

c ,A "=-°° J ° 

Here a number of comments are in order. First, in Eq.(^) is written for positive 
orbital momenta. This is sufficient because the negative ones can be obtained by 
changing the signs of a c and (3\ as can be seen from Eq.(||). Second, there are 
cancellations between contributions with different signs of the orbital momenta. 
For instance, there are even powers in 1/v in the uniform asymptotic expansion. 
They cancel!]. Also, contributions odd in f3\ cancel in the sum over the polar¬ 
izations. For large u, because of the definition Eq. (p7|) the difference of the two 
logarithms in Eq. (|57|) decreases as 1/iA In fact, under the sum over color and 
polarizations it decreases as l/v b . As a function of z it decreases as 1/z 4 . The 
integral over z decreases as 1/z/ 5 . So the convergence is quite well. For instance, 
it is sufficient to take a few orbital momenta to achieve a precision of several 
digits. 

In order to represent the results we introduce the notations 

£\v } z) = - V (in f^(iuz) - In fjf’ x)as (ivz)) (58) 

7T . ' ' 

c, A 

where the factor of 2 resulting from the two signs of the orbital momentum was 
taken into account and 

OO 

£ i {y) = J dz z £ { (is,z) (59) 

O 

3 This is taken into account already in the coefficients X n j and Z n j listed in the Appendix. 
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30 K 


Figure 1: The s-wave contribution, S^k), multiplied by k A . 


where the factor of 2 results from the two sign of the orbital momentum so that 


£' = £ (so) 

v=\ 


holds. 

In a similar way the contribution of the s-wave (v = 0) can be treated. Actu¬ 
ally, it is much easier because there is only one integration and no sum. We put 
m = 0 in Eq.([37]), integrate by parts and write it in the form 


OO 


with 



As an example let us consider 


£ 0 f = j dkk £l(k) 

0 

(61) 

yin 

(62) 

c, A 


the profile 


n[r) = tanh 2 (r). 

(63) 


It vanishes for r —► 0 together with its first derivative and it approaches its value 
at infinity exponentially fast. 

Fig. 1 shows the function S^k) multiplied by k A . It can be seen that its 
decrease is confirmed by the numerical results. 

In Fig. 2 the functions £ f (z/, z) multiplied by z 3 u A ' 5 are plotted showing that 
the expected behavior for large z is confirmed by the numerical evaluation. The 
factor z 3 is taken to demonstrate the fast asymptotic decrease and the factor z/ 4 5 
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Figure 2: The functions £ f (y, z) multiplied by z 3 z/ 4,5 for the first ten orbital 
momenta. 

in order to have the curves for all orbital momenta (u = 1,..., 10) to fit into one 
figure. 

In Fig. 3 the contributions £ f (z/) from the higher orbital momenta (v = 
2,..., 10) to the vacuum energy multiplied by u 5 are shown. The corresponding 
sum in Eq.(|6l|) is fast converging and in fact it is sufficient to take some few first 
orbital momenta only. As an illustration we show the corresponding numbers, 

£q=-0.0482214, 

v= 1 2 3 4, (64) 

u 2 S i (u)= -0.00587755 -0.00039075 -0.000156122 -0.0000753351. 

It should be remarked that the observed decrease of the subtracted Jost func¬ 
tions is a strong check confirming that the formulas for both, the Jost function 
itself and for its asymptotic expansions, are correct. This is because these for¬ 
mulas are obtained in completely different ways. The first one is a numerical 
solution of the differential equation and the second one results form the asymp¬ 
totic expansion which we obtained from the Lippmann-Schwinger equation. 

Further, it should be noticed that in each case the decrease is by one power 
better than it follows from what we subtracted. This is because each second 
power disappears under the sum of color and polarization. 

As already mentioned, it turns out that for v = 0 and for v — 1 there 
are bound states. They manifest themselves as poles of the logarithm of the 
Jost function and, as a consequence, of the functions £q (k) and £ f (z/, k). These 
functions are shown in Fig. 4. 
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Figure 3: The contributions S i (u) to the vacuum energy multiplied by u 5 from 
the higher orbital momenta 


4 The asymptotic part of the vacuum energy 


Next we have to calculate the asymptotic contributions to the vacuum energy. 
Actually, the problem is to remove the auxiliary mass m and perform the analytic 
continuation to s = 0. We start with £q s , Eq.(^) and insert there the asymptotic 
expansion (£f^). For s > | we perform the limit m —> 0. The contribution from 
In /q s ('</c) which is proportional to 1 /k vanishes. After that we can put s = 0 and 
are left with the finite integral 



dk k 2 


d 1 
dk (k 2 + l) 3 / 2 


-2 


and obtain 


, as = J_^[dr(m P 2 (r)' 
0 —8n x r \ 16r 3 

c,A o \ 


8r 3 


(65) 


where P(r) is given by Eq. ([40|) . The sum over color and polarization can be 
carried out using (j^) and (|I0|), 


El = 2 N, E«c = 2 N, 


c, A 


c, A 


= 4 (SU(2)), 

c, A 


c, A 

We obtain 


Zai = 9/2 (SU(3)), E » 2 C (^-Pl) 

- ' c ,A V y 


— 11N 
6 ‘ 


pas _ 

On — 


OO 

N J y (VO) - j^. 4 (r) - 8r 2 (/i'(r)) 2 ^ 


32tt 


( 66 ) 


( 67 ) 


for the s-wave contribution to the asymptotic part of the vacuum energy. 
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Figure 4: The functions £g(k) and £ f (l ,k) for small k. The poles indicate the 
positions of the bound states. 


The contributions from the higher orbital momenta require more work. We 
start from rewriting £*“, Eq.(^), using the Abel-Plana formula (see Appendix, 
Eq. (|98|)) in the form 


with 


£»as, a 
£as, 0 

£ aS ’ b 


c, A 


£ as = j2\ s 33 ’ a - -s 


as, 0 


+ £ 


as, b 


( 68 ) 


C, 

—47T 

Cs 
—4(/r 
C a 

— 4(7T 


r du r dk (k 2 -m 2 Y S 4-ln/i c,A)as (^), 
JO Jm V J ok 

roc / x 1—s f) 

lim / dk(k 2 -m 2 ) —In &’ X)as (ik), 

—o Jrn ^ > dk Jv K h 



1 


X- 

i 


( ln /(t)' ) °‘(*) - ( cx -)) 


(69) 

(70) 

(71) 


where (c.c.) denotes complex conjugate. It should be noticed that /^ c,A ) as (ifc) 
is real so that (c.c.) changes in fact the sign of iv only. In Eq.(^8|) a factor of 
2 from the negative orbital momenta was taken into account. The contribution 
from v — 0 deserves a remark. In fact, Eq.([J5[) is valid for a finite /(0) only. 
This is ensured in our case by the auxiliary mass m. So we have first to perform 
the limes v —> 0. After that the limit m —> 0 produces from £ as ’ 0 a divergent 
contribution proportional to 1/m which is, however, canceled by a corresponding 
contribution from £ as ’ b . 
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Now we use the uniform asymptotic expansion of the logarithm of the Jost 
function as given by Eq.([44|). So we insert it into £ as ’ a and interchange the orders 


of integration. For each contribution of the sum over n and j we use formula (IOC) 
given in the Appendix and obtain the representation 


£as, a _ 


1 + s(2 In 2 — 1) 


8n 


m 


2s 


dr 


r (2-s)r(s) / 


J 3j 


i — 1 ’ 


(72) 


0 c, A j =3 2 

where contributions proportional to powers of m are dropped. From Eqs. (0) 

-11JV 


and (^ 6 ]) we obtain 


EEt 


"3 j 


l _ 1 
c, A j =3 2 1 


6 


irn'{r)y 


In the limit s —> 0, £ as ’ a takes the form 


£*as, a _ 


— 11IV 
487T 2 


A (—f- 2 In 2 — 1 — 7 — 21nm + 0(s 


(73) 


(74) 


with A given by Eq. 

Next contribution to be considered is £ as ’ °, Eq.(f70|). We insert In from 

Eq. (|44|) and use formula (p9|) from the Appendix, 


3 3 n 


e “-° = fiSi/vNN £2* 

q c,A n= 1 j=3n 


X 


r(2 - s)r(s + \ - 1) m 2 ~ 2s ( v) 3 ~ n 

2\ s +2 _1 


(■ mr) j r(|) 


! + (—) 

\mr J 


(75) 


Here in the limes v —> 0 only contributions with j — u remain (note j > n in 
Eq. (pI2|) ). After that we consider the limits m —>■ 0 and s = 0. We obtain 


dr 




(76) 


c.A 


The contribution from n = 2 is odd under j3\ —> —/?a and does not contribute. 

It remains to consider £ as ’ b , Eq. ([H|) . Again, inserting ln/ as (ifc) from Eq.(fPI|) 
and using formula (|99| ) from the Appendix, 


£ 


as, b _ 


C a ) 


roo f) p fj,p 3 3 n 

L * V 

o n=l]=n 


An 


w 2—2s r(2 — s)r(s+ § — 1) 
xm -;-- £, n j(mr) 


{mr)i r(|) 


( 77 ) 
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with the notation 


E ” ;(X) / l-e^i /,,„J 


(t) + 1 


s+j/2-1 


— C.C. 


holds. In Yi n j(x) there is in fact no contribution from v < x because j — n is 
even. It can be rewritten in the form 


£'nj (*^) 


r°° dv (—l)( n 1 )/ 2 cos 7is v j 


1-e 2 


(s) 2 - 0 


s+j/2-l 


We have to perform the analytic continuation to s = 0. For n — j — 1 we simply 
can put s = 0 and obtain for small x (this is because we need m —*• 0 ) 

Sn ^ = 12x +0 ^' ( 80 ) 

For n — j — 3 we also can put s = 0 and obtain 


= x 


i- 00 dv 1 

e 1 — e 2 ™ \/z/ 2 — a: 2 


This integral is convergent. But for x —► 0 this form is not convenient. So we 
integrate two times by parts. We obtain first 


£ 33 ( 0 ;) =--- 2 / d/y arctan 

V ; 1 -e 2 ™ J x 


/ 2 — x 2 \ 1 — e 


and subsequently 


£.i.i(V) = 


1 -e 2 


+2 / d/y I /y arctan 


+ a; (-7T + 2 In 2x 2 ) (^ 2 ^) 
an —^== + x [in 2 x + In [/y + \/v 2 — x 2 ) ) ^ 


1 -e 2 


Here we can tend a; to zero under the sign of the integral and arrive at 

£33 (a) = \ + x ( lna; + S s) + 0(x 2 ) 


E --i +2 f* (tzW)' 


= -1.36437. 


For n — 3, j — 5 we proceed in a similar way. In order to perform the analytic 
continuation to s = 0 we integrate by parts and obtain 


Eas fx) = — 2ar 


r°° dv 


lx — x 


2 — x 2 V 1 — e 
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For x —» 0 we have to integrate by parts once more and get 

£35(20 — x 3 (l nx + £5) + 0 (x 4 ) 

with 


E 5 = 2 / dv \n(2v) 
Jo 


1 — e 


2nu 


= -0.364378. 


( 87 ) 


( 88 ) 


Finally we have n = 3, j = 7. For s = 0 we integrate two times by parts and 
obtain 


2 r c 

£37(20 = -3 X 5 J x 


dv 


/\ ' 


v 


\fv 2 — x 2 \v yl — e 2nu ) 
For i-^Owe have to integrate by parts once more and get 

T, 37 (x) = x 5 (In x + E 7 ) + 0(x 6 ) 


with 


£7 = 




/\ n 


= -0.121459. 


(89) 


(90) 


(91) 


Having now the appropriate representations for E n j we can insert them into 
£ as ’ b , Eq.([77|). We start from considering n — j — 1. Inserting (^) we obtain 


£»as, b(ll) _ 


N 


dr 


48(/r Jo r 2 




(92) 


Next we look for the first contribution in the rhs. of Eq.(|84]). It gives a con¬ 
tribution to £ as ’ b , Eq.([77]) which is proportional to 1/m and it cancels just the 
contribution of £ as> °, Eq.(^). Next we consider contributions from E 3j - contain¬ 
ing the logarithm ln(rm). What is proportional to lnm can be seen to cancel 
exactly the corresponding contribution from £ as ’ a , Eq.([74|). The remaining part 
is 

-11N r°° dr 


£»as, b(ln) _ 


lnr //(r) . 


(93) 


247T Jo r 

Finally we collect the contributions from E j and denote them £ as ’ b ^ 3 b (j = 
3, 5, 7). They read explicitely 


pas, b(3j) = J_ f°° & y- y- Z 3j y 

4:71 Jo r 3 ■— j j /2 — 1 


(94) 


with Z 3 j given by Eqs.(P7|). Performing the sum over color and polarization we 
obtain 

£as, b(33) _ f (— Hp(r)/i / (r) + 6r/i(r)p"(r) — 6r(//(r)) 2 

24:7r Jo r 2 v 

-y/A’'V( r )) 1 
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£as, b(35) = NE 5 r°°dr 


£as, b(37) 


24tt Jo 72 ( 26 / i (' r )/ i, (' r ) - 6r/x(r)/i /, (r) - 5r(//(r)) s 

+iV/i 3 (r)/i , (r)) , 

-5NT, 7 r°° dr 


lQn 


l ^ /‘('VM- 


( 95 ) 


In this way the complete asymptotic part £ as from the higher orbital momenta, 
Eq.(|68|), of the vacuum energy reads 


£as _ £ as ’ t*(ll) _j_ £as, b(ln) _j_ ^as, b(33) _j_ £as, b(35) _|_ £as, b(37) (96) 

where the individual parts are given by Eqs. d92j), (|93|) and 


H(r) 

v - 

= 0 

v - 

= 1 

tanlT r 

K 0 

= 0.82862 

Kl 

= 0.16771 

tanh 4 r 

K 0 

= 0.58388 

Ki 

= 0.16379 

00 

4-^ 

K 0 

= 0.46125 

K 1 

= 0.15486 

r 6 /(1 + r 6 ) 

K 0 

= 0.742808 

Kl 

= 0.28539 


Table 1: The bound state levels k in Eq. fll6|) . 


5 Numerical Results 

The actual calculation have been performed for four prohles /i(r), see for example 
Table 1, using Mathematica. For the solution of the differential equation (55) 
the build in function NDSolve had been used. The working precision was taken 
to be 25 digits and the pure computation time was about 40 hours on a standard 
PC. The good convergence of the integrals and of the sums was already discussed 
in Sect. 3. The numerical integrations had been done with the build in function 
which calculated about 200 points for each integral. In this way the numerical 
effort is not very demanding. 

In the examples considered here there are bound states in the spectrum of the 
fluctuations. They result in an imaginary part in the vacuum energy, see Eq.(PT|). 
In the subdivision of the vacuum energy according to Eq. (|33|) the imaginary part 
resides in the finite parts, £q and £ { . Thereby the first bound state appears in the 
s-wave contribution £ q, the second in the p-wave (u — 1). The binding energies 
are shown in Table 1. These bound states do not depend on the gauge group 
SU(2) or SU(3) because they come from the color eigenvalues a c = ±1 in Eq.(|9|) 
which are common to both groups. 

In Table 2 the numerical results for the asymptotic parts of the vacuum energy 
£ as , Eq.(|35D are shown. It is seen that £ as ’ b ^ 33 ^ is the numerically dominating 


21 







N 

£as, b(ll) 

£as, b(ln) 

£as, b(33) 

£as, b(35) 

£as, b(37) 

£as 

2 

0.00823767 

0.16906 

fi(r) = tanh 2 r 
0.507441 -0.138928 

0.0150082 

0.208057 

3 

0.0123565 

0.25359 

0.753783 

-0.204451 

0.0225122 

0.311155 

2 

0.00360497 

0.00354893 

/i(r) = tanh 4 r 
0.272918 -0.0634807 

0.00656787 

0.138266 

3 

0.00540745 

0.00532339 

0.405012 -0.0928895 

0.00985181 

0.2064 

2 

0.00213291 

-0.0326859 

/i(r) = tanh 8 r 
0.190348 -0.0386886 

0.00388594 

0.0825228 

3 

0.00319937 

-0.0490288 

0.282582 -0.0564628 

0.00582891 

0.122963 

2 

0.00534584 

0.0228004 

/r(r) = r 6 /(1 + r 6 ) 
0.561348 -0.0998484 

0.00973955 

0.310584 

3 

0.00801875 

0.0342006 

0.833918 

-0.145444 

0.0146093 

0.464541 


Table 2: The constituents of the asymptotic part £ as of the vacuum energy ac¬ 
cording to Eq. ([96]) . 


contribution. Within it the dominating part is (//(r)) 2 over weighting that of 
Therefore there seems to be no chance to change the sign of 8^' b( ' 33 ^ 
taking another profile function n{r). 

In Table 3 the numerical results for the vacuum energy 8*™^ given by Eq. (|33| ) 
are shown. The dominating part is the asymptotic contribution £ as , followed by 
The finite contributions 8 { and 8q are considerably smaller. 

We conclude the representation of the numerical results by Table 4 showing 
the quantities entering the calculation of the minimum of the complete energy, 
namely A, Eq.(^6|), B , Eq.(^8|), the radius r 0 | min , Eq.fl30|) and the energy in the 
minimum, £| min , Eq.([H]) (the latter two quantities are given in units of Aq CD ). 


6 Conclusions 


We calculated the vacuum energy for a color magnetic vortex in QCD. We found 
that the renormalized vacuum energy calculated at r 0 = 1 is positive. As 

discussed in Sec. 2 the sign of the complete energy in the minimum £| min , Eq.(|3T|), 
is not affected by this sign but only its deepness. For a positive 8^nt it is less 
deep. At once the core radius becomes larger. In this sense fi(r) = tanh 8 r is the 
’best’ one among the considered profile functions. 

In general, the sign of the vacuum energy is hard to predict. In the given case 
the vacuum energy which is positive results from bosonic fluctuations 


whereas in the magnetic background calculated in || [H|, [|] it is negative resulting 
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N 

c as 

C 0 

T 38 £ { 0 

£ f 

CrenPV 
^ quant 

2 

-0.352761 

= tanh 2 r 
0.560818 -0.0482214 

-0.00659258 

0.153243 

3 

-0.526636 

0.83779 -0.0894967 

0.00379912 

0.273678 

2 

-0.0848931 

fi(r) = tanh 4 r 
0.223159 -0.0889963 

0.000655135 

0.0499246 

3 

-0.126305 

0.332705 -0.12675 

0.0101892 

0.178836 

2 

-0.0424693 

fi(r) = tanh 8 r 
0.124992 -0.0649512 

0.00630883 

0.0238805 

3 

-0.0631556 

0.186119 -0.0903911 

0.0161882 

0.113711 

2 

-0.188802 

yu(r) = r 6 /(1 + r 6 
0.499386 -0.08936 

) 

0.0452825 

0.266506 

3 

-0.280762 

0.745303 -0.137446 

0.0877988 

0.504254 


Table 3: The constituents of the vacuum energy according to Eq. (|33|) . 


from fermionic ones. However, these results cannot be compared directly because 
of different normalization conditions used (massless vs. massive case). Neverthe¬ 
less, one may wonder to what extend this is a general rule and what the quark 
contributions might change for a center-of-group vortex. 

The topic of stability of the vortex is to a large extend beyond the scope of 
the present paper. So we only remark that the bound states and the imaginary 
part found here make the whole configuration unstable with respect to creation 
of particles. This is the same instability as known from the homogenous color 
magnetic field. 

The bound states are well known in a homogeneous color magnetic field. 
In an inhomogeneous one their existence can be predicted from some simple 
arguments. First note that the background we are concerned with is in fact 
Abelian and the contribution proportional to //(r) in Eq.(|8|) is an additional 
magnetic moment. From [|22| it is known that for a balance between the lowest 
Landau level and an attractive interaction between the magnetic moment and 
the magnetic background field there are zero modes, one for each flux quantum. 
In |22[ this had been established for an electron but it is obvious that this holds 
in our case too. Here the magnetic moment is larger due to the spin of the gluon. 
The zero modes become more tightly bounded and turn into bound states. This is 
independent from the specific form of the background and it is what we observed 
for our profile functions The bound states found here are also in line with 

those in [|23| where an electron with gyromagnetic factor larger than 2 had been 
considered. 
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N 

A 

B 

n)| min A QC D 

^Imm/A-QCD 



H(r) = 

= tanh 2 r 


2 

0.864328 

0.148238 

11.1221 

-0.000324477 

3 

0.864328 

0.103024 

16.0032 

-0.000235092 



H(r) = 

= tanh 4 r 


2 

0.488359 

0.332642 

4.95644 

-0.000923169 

3 

0.488359 

0.072185 

22.8402 

-0.0000652096 



H(r) = 

OO 

cd 

+3 


2 

0.353234 

0.482925 

3.41403 

-0.00140737 

3 

0.353234 

0.0991947 

16.6211 

-0.0000890673 



/i{r) = r 

,8/(1 +r 6) 


2 

1.07484 

0.0692783 

23.7985 

-0.0000881308 

3 

1.07484 

0.0344776 

47.8201 

-0.0000327415 


Table 4: Summary for the energy of the vortex 


The calculations had been performed in an efficient way, using Jost functions 
and separating an asymptotic part thus leaving fast convergent expressions. The 
main numerical effort is for the finite parts. Here one has to perform the integra¬ 
tion over the radial momentum and the sum over the orbital momenta which are 
fast converging dne to the subtraction of the corresponding first few terms of the 
asymptotic expansions. For instance, in Eq. (p7|) we subtract up to O which 
is sufficient for convergence, ft turns out that the next power, \ is absent and we 
are left with an -^-behavior making the sum in Eq. (|57| ) fast convergent. The ab¬ 
sence of the -4-contribution is due to cancellations between the orbital momenta 
of different sign]]. Furthermore, this cancellation is related to the smoothness 
properties of the background. If the background has a jump, as, for instance, 
considered in || [24], [HJ], a contribution proportional to is present slowing 
down the convergence. Of course, one may think about over subtractions, i.e., 
including more terms in the asymptotic part of the Jost function, Eq.([4^). But 
there one needs more orders in the iteration of the Lippmann-Schwinger equation 
than done in [,6j. This had been done only for the particular example of the pure 
magnetic background in 
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In Sec. 5 it was found that the numerically dominating contribution results 
from £ as . It is of some interest to discuss the question whether £ as can be at¬ 
tributed to a orbital momentum. We start with the remark that the initial, 
regularized vacuum energy £ quant , (0), which appears after separation of vari- 

4 We note that this had been taken into account in the coefficients X nj and Z n j which are 
listed in the Appendix. 
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ables, is a sum over orbital momenta. This remains true also after the subdivision 
into finite and asymptotic parts in Sec. 3 and for the finite parts themselves. Also 
£p s can be attributed to the s-wave. ffowever, cannot. Here the reason is that 
in order to perform the analytic continuation to s = 0 we used the Abel-Plana 
formula and performed the integration over v after what it is impossible to iden¬ 
tify contributions from some v. In this way we cannot say from which orbital 
momentum the numerically dominating contribution originates. To some extend 
this discussion can be continued with the remark that it is the counter term who 
cannot be represented as a sum over orbital momenta. 
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Appendix 

The coefficients X n3 appearing in Eq.( 
a c 2 n(r) 2 


read 


*i,i = 

Ad,2 = 

*1,3 = 

* 2,2 = 

*2,3 = 

*2,4 = 


o, 

- (a c 2 n(r) J 


a 2 /i(r) 2 a 2 r /i{r) n'(r) 
2 2 : 
- (a 2 f3 x rn(r)n'(r)') , 

-5 ay 2 /i (r) 2 


+ a 2 r n(r) 


* 2,5 = 0 , 

* 2 ,6 = 2 a 2 n(r ) 2 , 


*3,3 = 


13 a c 2 /i(r) a c 4 /i(A 


16 


a 2 r /i(r) //(r) + 


a c 2 r 2 n'(ry 


a 2 P\ 2 r 2 n'(r) a 2 r 2 n(r) /i"(r) 

2 + 4 ’ 

*3,4 = -2>a 2 (3 x r n(r) n\r) + a 2 (3 x r 2 n'(rf + a 2 (3 x r 2 n{r) n"{r), 


*3,5 = 


—153 a 2 /r(r) 2 3 a c 4 /i(r) 4 27 a 2 r n(r) /i'( r ) 

16 + 4 + 4 
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5 a 2 r 2 n'{r) 3 a c 2 r 2 /j,(r)/x"(r) 

8 4 ’ 

X 3i6 = 4 a c 2 f3 x r /x(r) /d'(r), 

315 a c 2 /i(r) 2 5 a c 4 /r(r) 4 25 ay 2 r fi(r) n'(r) 


X 3.7 — 


16 


X;i.& — 0, 

—175 a c 2 ji(r ) 2 


^3,9 = 


16 


The non vanishing coefficients Z n] appearing in Eq.(44) read 

^1,1 = Oi 2 r n(r) n'(r), 

Z 2 ,3 = - (a 2 ^r/i(r)/r'(r)) , 


^3,3 — 


(a 4 r/r(r)V(r)) a 2 (-l + 2/? 2 ) r 2 n'(rf 


+ 


6 

a 2 r/i(r) ( —11/i'(r) + 6 r/i"(r)) 
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2 « \ 2 , 2 


Z 3j4 = a 2 /3r 2 ii'(r) + cr/3r/i(r) (—3//(r) + r//(r) ), 

a 4 r/i(r) 3 /i'(r) 5 a 2 r 2 / u , (r) 2 a 2 r/r(r) (—13/i'(r) + 3 r//'(r) ) 


■^3,5 — 


Z 3i6 = 4a 2 /3r/i(r)/i'(r), 
-25 a 2 r /r(r) //(r) 


^3,7 — 


(97) 


The Abel-Plana formula used in Sect.2 reads 


£/m=/<wm-±/(o)+ / ' *' 


i/=i 


1 - e 


2ttu 


(98) 


The following formulas are used in the text, 


r °°dk(k 2 -m 2 ) 1 ~ s ^ = -m 2 - 2sT ^ s)r(s+ 2 X ) (-) 


j-n 


dk 


r(f) 


1 + ( 


V 

mr 


2\ 


S+4-1 


(99) 


and 


poo /*oo 


du / dk(k — m ) 


2 2\l-s ^ V 


m 


2-2, r(2 - s)r(±±P)T(s + S=3) 


<9fc z/” 


2 (rm) 


n—1 


r(i/2) 


(100) 

with f = l/^/l + ( kr/v ) 2 . They can be easily derived, see also (C3) and (C2) in 
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